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Abstract 

■ We apply the theorems from the theory of Toeplitz determinants to calculate the 
asymptotics of the correlators in the XY spin chain in the transverse magnetic field. The 
asymptotics of the correlators for the XX spin chain in the magnetic field are obtained. 

<D ■ 

C/3 ; 1. Introduction 

m ■ 

7— i ■ The calculation of correlators in the one- dimensional exactly solvable models remains an 

interesting and important problem. An important exact results on the correlators for the XY 
quantum spin chain and the Ising model were obtained [l]-[3]. The goal of the present letter 

^ ■ is twofold. First, we present the calculations of the correlators in the limit of the XX spin 
chain in a non-zero magnetic field which was not considered previously. In particular, the 
asymptotics of the exponential correlator for the XX spin chain is found. Second, we use the 
generalized Fisher-Hartwig conjecture (for example, see [4]) and the method to calculate the 
asymptotics of Toeplitz determinants [5] to derive the known results [3] for the correlators in 
P-r the XY spin chain in a non-zero magnetic field in a simple way. 

■ In the present letter to calculate the asymptotics of the determinants we use the conjecture 
from the theory of Toeplitz matrices, [4], the generalized Fisher-Hartwig conjecture, which is 
based on the proofs of the original Fisher-Hartwig conjecture [6] including the constant in 
front of the asymptotics, in a number of the particular cases [7], [8] (see also [4] and references 

^ ■ therein). One finds the representation of the generating function f(x) of the Toeplitz matrix 
= M(i - j) = J^(dx/27r)e i ^ i - j '> x f(x) of the following form: 

f(x) = f (x)]Je ibr{x - Xr -^ n{x - Xr) \2 - 2cos(x - x r )) a \ (1) 

r 

where x G (0, 2n) is implied, the discontinuities (jumps and zeroes or the power-law singulari- 
ties) are at finite number of the points x r , and fo(x) is the smooth non- vanishing function with 
the continuously denned argument at the interval (0, 2ir). The function (1) is characterized by 
the parameters a r , b r at each point of discontinuity x r . In general there are several representa- 
tions of the form (1) for a given functions f(x). To obtain the asymptotics of the determinant 
one should take the sum over the representations (1) corresponding to the minimal exponent 

D(N) = e loN N^ a r- b r)E (2) 

Repr. 
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where E is the constant independent of N, 



E = exp £ fcZ fc Z_* n(/+(^))" ar+6r (/-(^))"° r " br ( 3 ) 

\fe=l / r 

TT i(x -x )\-K+b r )(a s -6 s ) „ G(l + <2 r + 6 r )G(l + CL r - b r ) 

M l S ^ V G(l + 2a r ) ' 

4 = f* n (dx/2ir)e lkx ln(f (x)), and the functions f±(x) are given by the equations 



—ikx 



fc>0 fc>0 

where G is the Barnes G-function [9], G(z + 1) = G(z)F(z), G(l) = 1. 

For the piecewise continuous generation function f(x) with the finite number of discon- 
tinuities which is not equal to zero (a r = 0) the equation (2) can be represented in a more 
convenient form [7]. Suppose the function f(x) has the finite number of discontinuities at the 
points x r , 

A r = ^(ln/(x r + 0)-ln/(a; r -0)), 

and assume that the function f(x) has the continuously defined argument and not equal to 
zero at the interval (— tt, it). Then the asymptotic of the determinant is given by the formula 
[7] 

D(N)= £ e l ° N N^ r *)E, E = exp(f: (kl k U - \ ^ A) H~g(\ r ), (4) 

Repr. \k=l V K r / / r 

where we denote l k = j^(dx/2^)e lkx hi{f(x)) and the function g(X) equals 

oo / \2\ 

^(A) = e^ 2 n(l + ^)e- A2/fc , (5) 



k=l 



where 7 = 0.577.. is Euler's constant. The constant E (5) in eq.(4) is a consequence of the 
formula for the Barnes G- function: 

G(l + z) = (2vr)^ 2 exp {-^^- - 7y) ft (l + $ exp(-^ + z^/2k), 

in such a way that g(X r ) = G(l + z)G(l — z) for A r = iz. 

In ref.[7] the theorem was proved for an arbitrary number of the discontinuities of the 
imaginary part of the magnitude less than 1/2, |A r | < 1/2, however, there are many reasons 
to believe it to be true also in the case |A r | = 1/2 [4]. For the case of an arbitrary single Fisher- 
Hartwig singularity the conjecture (2) was recently proved in ref.[8]. See ref. [4] for the complete 
list of the cases for which the rigorous proof of the conjecture (2) is available. Schematically 
the rigorous proofs in the particular cases [7] [8] go as follows. Suppose that equation (4) 
is fulfilled for some functions fi(x) and fi{x) of the class of piecewise continuous functions 
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with continuously denned argument. Then the equation (4) is fulfilled for the function f(x) = 
fi(x)f 2 (x). Thus it is sufficient to prove eq.(4) for the smooth function with the continuously 
defined argument, in which case it is reduces to the well known strong Szego theorem and for 
the singular function of the form (1) f(x) = (1 — z) a (l — 1/z) 13 , z = e lx , in which case the 
asymptotics is known exactly. 

2. General form of the correlators for XY spin chain. 

Let us review the known calculations [l]-[3] of the equal-time correlators of the operators 
S^ 1 for the XY spin chain at non-zero transverse magnetic field. The Hamiltonian of the XY 
spin chain has the form 

H = - £ ((1 + 7)^1 + (1 - l)SfSf +1 - hSf) , (6) 
i=i 

where S a = ^a a , a = x, y, z, < 7 < 1 and the periodic boundary conditions are implied. 
Performing the well-known Jordan- Wigner transformation (S^ = S x ± iS y ): 

S H- = e ^(*) a + = exp(5:n I )a+ (7) 

Kx 

we obtain the following Hamiltonian written down in terms of the fermionic operators: 

L L 
^ i=l i=l 

where the boundary terms which are not important in the thermodynamic limit are omitted. 
The Hamiltonian is diagonalized with the help of the following Fourier transform in the sectors 
with an odd number of particles: 

a t = -7rY, eikXa t, fc = -7^, -7T<k<7i. 
V L x L 

The Hamiltonian (8) is easily diagonalized in the momentum space. Introducing the column 
V'fc" = ( a ti a -k) for k G (0, 7r) the Hamiltonian is represented as 

// E'v\//,'.v, M k =( 6k tAk ), 4> k =( a *), k>0, 

where the notations e k = cos(k) — h, = •ysm(k) are used and the additional constant 
term e k is implied. It is convenient to redefine the operator at k — > id-k- ^is case the 
Hamiltonian is diagonalized by means of the canonical transformation 

a t = c k&ik + Sfe«2fc, a-k = c k a 2k - s k af k , (9) 



H = --J2 ( a t a i+i + l a t a t+i + h - c ) + hY.attti, (8) 
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where k > 0, c 2 . + s\ — 1, c k = cos(0fc), Sfc = sin(0fc), and the angle <p k and the Hamiltonian 



are: 



tg(20 fe ) = H=-J2 (at k a lk - a\ k a 2k + l) , E k = ^/(cos(fc) - hf + 7 2 sin 2 (A;). 

6fc fc>0 

Consider the equal-time correlators for the XY spin chain: 

G(x) = (0\S+JSr\o). G(x) = (0\S+JS?\0). 

Using the Jordan- Wigner transformation S£ = exp(i7iJ2i<x n i) a ti the correlation functions 
G(x), G(x) can be represented as the following averages over the ground state of the Hamil- 
tonian (8), 

G(x) = (0\ate inN{x) a \0), G(x) = -(0|a+e^ (:c) a+|0), (10) 
where N(x) = Z)f=i n i- Introducing the operators, anticommuting at different sites, 

Ai = at + a u Bi = af-ai, A i B i = e i ™\ 

where rij = afdi - is the fermion occupation number, with the following correlators with 
respect to the vacuum, 

(BiAj) = BA{i - j), (AiAj) = 0, (B.B,) = 

one obtains the following expression for the bosonic correlator: 

1 X — 1 1 X — 1 

G(x) = -(0\B X n(AA)A)|0) + -{0\BoYl(ABi)Ao\0) = (S y x S y ) + (5£S*>, (11) 
4 i=i 4 i=i 

where the terms (S^.Sq), (S*Sq) correspond respectively to the two terms at the left-hand 
side of this equation. At zero magnetic field each term in this expression is represented as a 
product of the two Toeplitz determinants [1]. The asymptotic behaviour of G(x) is determined 
by the form of the average BA(x). To calculate this function one should use the expressions 
for the averages of two operators af, which follow from the canonical transformation (9): 

(a + an \ = r<0L ( e ^ + e -**) ( I + 1 cos(k)-h \ 

Jo 2ttV Jy 2 2/ {cos{k) _ h) 2 + 1 2 sin 2 {k) )> 



sin(/c) 



2 v /(cos(A;) -/i) 2 + 7 2 sin 2 (A;) / 



where the sums have been replaced by the integrals in the continuum limit. From the expres- 
sions (12) the following expression can be obtained: 



which gives the following generation functions fi(x), f2(x) for the correlators (S%Sq), (S^Sq) 
respectively for the XY spin chain in the magnetic field: 

r i \ Tix- i \ Tix fcos(k) - h + rysin(£;)\ 1/2 

f l2 ( x ) = ensign a; = e Tlx }-f , , . 13 

' \cos(k) — h — -ry sm(fc)/ 

Finally, let us obtain the determinant representation for the exponential correlator 

G a (x) = (e iaN ^), (14) 
where a is an arbitrary parameter and N(x) = J2i=i n i, introducing the operators 

Ai = e ia/2 a+ + a h B t = a+ + e ia/2 a h A { Bi = e ian \ 

for an arbitrary a we represent the correlator (14) as the average (rii=i( a i-^))- Since for a ^ ir 
the averages (AjA,) and {BiBj) are not equal to zero at 7 7^ 0, we obtain the correlator (14) 
as a Pfaffian, which can be expressed through the Toeplitz determinants. For the case a = n 
we have (AiAj) = (BiBj) = and the average can be represented as the Toeplitz determinant 
det^ ({AiBj)), i,j — 1, . . .x, At 7 = the generating function for this determinant equals 
f(x) = sign(|x| -p F )- 

For the practical calculations one should represent the functions (13) in the following form: 




A, 2 (-)= --f— --f— , (15) 



where the factor e l2x in the parenthesis corresponds to the correlator {S^Sq) (to the function 
f2(x)) and the parameters Ai 5 2 are given by the expression: 



Ai,2 = Y^{ hT V^ 2 + 7 2 - l) • (16) 

The behaviour of the correlators is determined by the values of the parameters Ai^, which 
obey the relation A1A2 = o, where the parameter a < 1 equals a = (1 — 7)/(l + 7). In 
general, for 7 > one should distinguish four different cases. 1) h 2 + 7 2 < 1. In this case 
the parameters A ii2 = \fae ±%<i ' are the complex numbers. 2) h 2 + 7 2 > 1 but h < 1. In this 
case the parameters Ai ; 2 are real and obey Ai < A2 < 1. 3) The line h — 1 at arbitrary 7: 
Ai = a < 1, A 2 = 1. 4) The region h > 1 where the real parameters Ai j2 are Ai < 1 < A 2 . The 
asymptotic behaviour of the correlators is different in this four regions. 



3. Correlators for the XX spin chain. 

In this section we calculate the correlators for the isotropic XX spin chain (7 = 0) at non- 
zero magnetic field. For this purpose it is sufficient to use the equation (4). Let us present 
the generating functions f(x) for the initial determinants of the x x x matrices obtained for 
different correlators. We construct the functions with the continuously defined argument at 
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the interval (— ir, ir) which is important for the application of the equation (4) (see [7]). One 
can use the following functions. First for the correlator (14) at a < vr one should use the 
function 

/(x) = (l,e fa ,l) (17) 

where the three number in parenthesis denote the values of the function at the intervals 
(— 7r, — pf), (—Pf,Pf) and QoF,vr) respectively and the unity in eq.(17) denotes the phase 
equal to zero, the function 

f(x) = (1,-1,1) (18) 
for the same correlator at a = vr and the functions of the form 

f 1 (x)=e^ x (e^,l,en, / 2 (x) = e^e**, 1, e"*), (19) 

with the continuously defined argument for the two terms in eq.(ll) respectively. Let us note 
that for the case a > vr in eq.(14) the function can be chosen in the form (e l27r , e ta , e t27T ), which 
shows that the correlator (14) is in fact the 2vr - periodic function of the real parameter a. 

The function f(x) (17) has the continuously defined argument and the two jumps of the 
size ia and — ia at the points X\ = —vr/2 and x 2 = vr/2. Substituting this functions into the 
equation (4) we obtain at the half-filling the asymptotics: 

G a (x) = (e^)> = e^' 2 ^^ (,(A)) 2 , (20) 

where A = a/2n. Evaluating the product in g(X) for the function given by the equation (5) 
we obtain the expression: 

i / / ^ \ \ x2 rdtf t 1 (sh(At/2) 2 \ x2 a 2 

For the correlator G a (x) at a = rr one should use the conjecture (2). The function f(x) = 
(1, — 1, 1) (18) allows for two different representations in the form (1) with two equal exponents 
~ l/x 1 / 2 : 1) bi = 1/2, b 2 = -1/2, f (x) = e™l 2 and b x = -1/2, b 2 = 1/2, f (x) = e"^ 2 , with 
X\ = —vr/2, x 2 = it/2. Thus for the asymptotics one should take the sum of the two terms 
which gives the following result: 

G n (x) = (e ipFX + e ~ ipFX ) G Q (x) = cos{p F x)2G {x), (22) 

where pp = vr/2 and 

Go{x) = ^-L(,(l/2)) 2 . 

The correlator G n (x) equal to zero for x - odd and to ±2Go(x) for x- even in accordance with 
the general properties of the average G n (x) = (Y\ X =1 {1 — 2raj)) (real function). The same result 
for Gn can be obtained using the reduction of the determinant to the square of the half-size 
determinant. At a ^ vr the first term e ippX in eq.(22) corresponds to the expression (20) while 
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the second term gives the subleading asymptotics. One can also obtain the asymptotics of the 
correlator (22) at the arbitrary p F ^ tt/2: 

G n (x) = (e*** + e-^)-^ T =i_i=(^(l/2)) 2 ! (23) 

x Jsm(p F ) V z 

which also have two different exponents and at arbitrary p F is the real function of x. Note 
that the same results could be obtained from equation (4) by taking the sum of the two terms 
corresponding to the functions (1, e l7r , 1), and (1, e _Mr , 1). 

The spin-spin correlator (10) using the functions (19) for the two terms we obtain from 
eq.(4) the asymptotics 

= ^^GKl/2)) 2 , (24) 
where the constant <?(l/2) = y/ir(G(l/2)) 2 , can be represented in the form: 

1*07(1/2)) = 




(ch(f)) s 



Note that for the functions (19) there is only one representation (1), b\ = b 2 = 1/2 and 
h = b 2 = —1/2 respectively so the conjecture (2) leads to the same results. 

Finally, one can use equation (11) and the functions (19) and the discontinuities at the 
points xi t 2 = ±Pf to obtain the asymptotics of the correlator G{x) for the XX spin chain in 
the magnetic field in the form: 



G{x) = ^^-|^(,(l/2)) 2 , (25) 

which differs from the well known asymptotics (24) by the simple factor (sm(p F )) 1//2 = (1 — 
h 2 ) 1 / 4 . Note that both eq.(2) and (4) predict the real positive function G(x) without the 
oscillating factors in agreement with the observation that the ground-state wave function of 
the Hamiltonian (6) at 7 = is the positive function of the coordinates. 

4. Correlators for the XY spin chain. 

Let us start with the calculations of the asymptotics of the ^-correlator G xx (x) = (S x Sq) 
for the different values of the parameters 7, h (Ai j2 ). In Section 2 this correlator was repre- 
sented by the determinant corresponding to the function fi(x) (15). At h < 1 we have for the 
parameters Ai^ the condition | A.x,2 1 < 1, and this function is the smooth function with zero 
winding number, so the strong Szego theorem can be applied. Thus in the whole region h < 1 
the correlator has the long-range order which can be easily calculated: 

G xx (x) = \(1 - a 2 ) 1 / 2 (l - h 2 fl\ a = ^—L, h<\. 
v ' 4 V ' v ' I + 7 
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Next, let us consider the correlator at the line h — 1. At this value of h the parameter A2 = 1 
and the function fi(x) takes the form: 

One can see that this function has the form (1) with the function fo(x) given by the first 
factor in the last equation and the single jump discontinuity at z — 1 with the parameter 
b = —1/2. The calculations according to the formulas (2), (3) yield the result 

1 -Y 3 / 4 

GXX W = ^^Ti } ^ (G(1/2))2 ' h = x > 

in agreement with the result of ref. [3]. Let us turn to the case h > 1. In this case A2 > 1 and 
the function takes the following form: 

To calculate the determinant we use the method proposed in ref. [5], which allows one to 
use the general Fisher-Hartwig conjecture (2), (3). Representing the matrix elements as an 
integrals in the complex plane over the unit circle, we deform the contour of integration to the 
circle \z\ = 1/A which is the outer boundary of the annulus A < \z\ < 1/A. This procedure is 
correct since the function (26) is an analytic function inside this annulus and reduces to the 
substitution z — > z/X in eq.(26). As a result we obtain the function of the form (1): 

1 = 1 , _3 

(l-xyzy/^ a 4' 4 

with the singularity at the point z — 1. Applying the conjecture (2) we obtain the result for 
the case h > 1: 



G xx (x) = 




h>l, 



in agreement with the corresponding formula in ref. [3]. 

Let us proceed with the calculations of the ^-correlator G yy (x) = (S^Sq). The behaviour 
of this correlator in the four regions in the space of parameters is completely different. We 
begin with the basic region in the space of parameters where both 7 and h are sufficiently 
small: ti 2 + •y 2 < 1. In this case the parameters Ai ; 2 are inside the unit circle at the complex 
plane: 




In order to obtain an appropriate representation (1) we deform the contour of integration from 
the unit circle to the circle of the radius y/a, which amounts to the substitution z — * \faz in the 



Hz) = x 



1 - XiX/z 
1 - (VAk 
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function f2(z) (15). One can see that the resulting function f(z) has the two singularities at 
the points z = e ±l( ^. One can prove that the only possible way to represent this function in the 
form (1) is as follows. Namely, we single out the two equivalent Fisher-Hartwig singularities 
at the points x = ±0. Thus we have the representation: 

f (z) = a(l - e^azY x l\\ - e^az)" 1 ' 2 , a x = a 2 = \, h = b 2 = ^ 

Substituting this function into the equations (2), (3) we obtain the correlator: 

G m (x) = a x - f^sin(0)(l-a)' 1 / 2 (l + a 2 -2acos(20))- 1 / 4 ") . 

Next, let us consider the case h 2 + 7 2 > 1, but h < 1, when the parameters Ai and A2 are real 
and both are inside the unit circle, Ai < A2 < 1. Taking the contour to be the inner boundary 
of the region A 2 < \z\ < 1/A 2 , where the function f2(z) (15) is an analytic function, we obtain 
the representation: 

k[z) ~ X2 \l- (X^z ) (l-A|z)V2' a ~l> 4' 
which gives the following correlation function: 

G^{x) = -1 (XIY -1(1 - A 2 ) 1/4 (l - A 2 )- 3 / 4 (l - A 1 A 2 )-V 2 (i - Ai/A 2 ) _1 , h < 1. 

X° y ' 07T 

This expression is identical with that given in ref . [3] which can be easily established using the 
following useful relations: l-AiA 2 = l-a = 27/(1+7) and (l-A^l-A 2 ,) = 4(l-/i 2 )/(l+ 7 ) 2 . 

As in the case of the xx-correlator, at the line h — 1 no special deformation of the 
integration contour is required. The parameters Ai j2 are Ai = a, A 2 = 1 and from the function 
/ 2 (x) (15) one can easily obtain the representation: 

f (z)= '-) , a = 0, b = -, 

\l-az J 2 

Calculating the constant (3) we get: 

G™{x) = -^^7- 5/4 (l + 7)(G(l/2)) 2 . 

Finally, let us consider the region h > 1 (A 2 > 1). In this case the calculations are analogous 
to the corresponding case for the xx-correlator. Introducing the parameter A = 1/A 2 < 1, we 
observe that the function is an analytic function in the annulus A < \z\ < 1/A. Deformation of 
the integration contour to the inner boundary of the annulus, z — > Xz, gives the representation 
(1) in the form: 

1/2 1 . 5 



4 4 
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Applying the conjecture (2) we obtain the result: 



G yy (x) = 




in agreement with the corresponding formula in ref.[3]. Let us note that in all the cases 
considered in this section (except the yy- correlator at h 2 + 7 2 < 1) we have obtained the 
functions with the single Fisher-Hartwig singularity. In the case of the single singularity the 
asymptotics (2), (3) was proved in ref.[8]. In conclusion, it is easy to verify that in all the cases 
considered, there is only one possible way to choose the contour of integration (inner or outer 
boundary of the annulus). In fact, in all the cases at the opposite boundary the condition 
a r ± b r G Z~ is fulfilled. Formally, in this case we have the constant E — 0, which suggests, 
that these special cases should be excluded from the general statement of the conjecture (2), 
(3) since in this case the matrix elements corresponding to the pure Fisher-Hartwig singularity 
does not exist. 
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